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Abstract 


The  acoustic  prediction  model  called  Bellhop  continues  to  be  enhanced  to  more  closely  fit  the 
requirements  of  DRDC  Atlantic’s  Environment  Modeling  Manager  (EMM).  This  version  4 
contains  both  passive  and  active  algorithms.  In  this  version,  linear  range  interpolation  of  the  SSP 
and  curvilinear  interpolation  of  the  bathymetry  are  added  as  input  choices,  and  an  additional 
output  of  the  sampled  SSP  is  provided.  The  major  differences  between  BellhopDRDC  and  the 
web  version  dated  May  2008  lie  in  the  input  data  and  file  formats  that  have  been  altered  to  satisfy 
the  requirements  of  the  controlling  programs  within  the  Environment  Modeling  Manager.  This 
document  provides  a  users  guide  to  the  running  of  the  active  and  passive  versions  of  the 
BellhopDRDC_v4  program  and  the  boundary  loss  program,  and  describes  some  plotting  routines 
available  for  viewing  the  prediction  results. 


Resume 


Les  ameliorations  du  modele  de  prevision  acoustique  Bellhop  se  poursuivent  afin  de  mieux 
P  adapter  aux  exigences  du  progiciel  de  gestion  de  la  modelisation  de  l’environnement 
Environment  Modeling  Manager  (EMM)  de  RDDC  Atlantique.  La  version  4  contient  des 
algorithmes  passifs  et  actifs.  Dans  cette  version,  V interpolation  lineaire  de  distance  du  profil  de 
vitesse  du  son  (PVS)  et  1’ interpolation  curvilineaire  de  la  bathymetrie  sont  ajoutes  aux  choix 
d’ intrants  et  une  capacite  de  sortie  supplementaire  sur  le  PVS  echantillonne  est  offerte.  Les 
principales  differences  entre  la  version  RDDC  du  Bellhop  et  la  version  sur  le  Web  en  date  de  mai 
2008  resident  dans  les  formats  des  donnees  d’ entree  et  des  fichiers,  qui  ont  ete  modifies  pour 
satisfaire  aux  exigences  des  programmes  de  commande  de  1’ Environment  Modeling  Manager.  Le 
present  document  constitue  un  guide  d’utilisation  des  versions  active  et  passive  de  la  version  4 
RDDC  du  Bellhop  ainsi  que  du  programme  de  perte  de  transmission  aux  limites,  et  decrit 
certaines  routines  de  tra9age  permettant  de  visualiser  les  resultats  des  previsions. 
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Executive  summary 


Users  Guide  to  BellhopDRDC_V4:  Active  and  Passive  versions 

McCammon,  D.F. ;  DRDC  Atlantic  CR  2010-134;  Defence  R&D  Canada  - 
Atlantic;  October  2010. 

Introduction 

Bellhop  is  a  computer  program  created  by  Dr.  Michael  Porter  that  computes  acoustic  fields  in 
oceanic  environments  via  Gaussian  beam  tracing.  The  environment  consists  of  an  ocean  that  may 
have  range  variations  in  the  sound  speed  profile,  the  bottom  loss,  and  the  bathymetry.  Two 
programs  were  created  for  use  with  the  Environment  Modeling  Manager  (EMM),  a  tactical 
oceanography  tool  for  naval  planning  and  operations.  The  first  is  the  passive  acoustic  version 
named  BellhopDRDC_ray_TL_v4.  The  outputs  from  this  program  include  transmission  loss 
(coherent,  semi-coherent,  and  incoherent)  and  ray  traces.  The  second  is  the  active  version  named 
BellhopDRDC_active_v4,  which  outputs  the  arrival  tables,  reverberation  time  series,  target  echo 
time  series,  and  the  signal  excess  versus  range. 


Results 

The  major  differences  between  BellhopDRDC  and  the  web  version  dated  May  2008  lie  in  the 
input  data  and  file  formats  that  have  been  altered  to  satisfy  the  requirements  of  the  controlling 
programs  within  the  Environment  Modeling  Manager.  This  document  provides  a  User’s  Guide  to 
the  running  of  both  the  active  and  passive  versions  of  the  BellhopDRDC_v4  program,  and 
describes  some  plotting  routines  available  for  viewing  the  prediction  results. 

Significance 

The  Environment  Modeling  Manager  is  a  sophisticated  tactical  oceanography  system  being 
developed  to  aid  naval  planning  and  operations.  It  provides  tactical  decision  aids  with  accurate 
and  consistent  predictions  of  acoustic  conditions  and  target  detectability.  The  Bellhop  software 
package  is  the  heart  of  the  system  that  provides  the  necessary  acoustic  predictions  to  client 
programs. 

Future  plans 

It  is  intended  to  continue  enhancing  the  Bellhop  program  with  more  accurate  models  that  are 
effective  from  an  operational  standpoint. 
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Introduction 

Cree  par  M.  Michael  Porter,  le  programme  informatique  Bellhop  permet  de  calculer  les  champs 
acoustiques  en  milieu  oceanique  par  traqage  de  faisceaux  gaussiens.  Le  milieu  est  un  ocean  qui 
peut  presenter  des  variations  de  la  portee  acoustique  dans  le  profil  de  vitesse  du  son,  de  la  perte  au 
fond  et  de  la  bathymetrie.  Deux  programmes  ont  ete  crees  et  seront  utilises  avec  V Environment 
Modeling  Manager  (EMM),  un  outil  d’oceanographie  tactique  pour  la  planification  et  les 
operations  navales.  Le  premier  est  une  version  acoustique  passive  appelee 
BellhopDRDC_ray_TL_v4.  Les  resultats  de  ce  programme  comprennent  l'affaiblissement 
acoustique  (base  sur  la  sommation  coherente,  semi-coherente  ou  incoherente)  et  les  traces  des 
rayons.  Le  deuxieme  est  une  version  active  appelee  BellhopDRDC_active_v4,  qui  permet 
d’obtenir  les  tables  d’arrivee,  les  series  de  temps  de  reverberation,  les  series  de  temps  d’echo  de 
cible  et  l’exces  de  signaux  par  rapport  a  la  portee. 


Resultats 

Les  principales  differences  entre  la  version  RDDC  du  Bellhop  et  la  version  sur  le  Web  en  date  de 
mai  2008  resident  dans  les  formats  des  donnees  d’entree  et  des  fichiers,  qui  ont  ete  modifies  pour 
satisfaire  aux  exigences  des  programmes  de  commande  de  V Environment  Modeling  Manager.  Le 
present  document  constitue  un  guide  d’utilisation  des  versions  active  et  passive  de  la  version  4 
RDDC  du  Bellhop  et  decrit  certaines  routines  de  traqage  permettant  de  visualiser  les  resultats  des 
previsions. 

Importance 

Id  Environment  Modeling  Manager  est  un  systeme  perfectionne  dans  le  domaine  de 
l’oceanographie  tactique  pour  le  soutien  a  la  planification  et  aux  operations  navales.  II  represente 
un  outil  d’aide  a  la  decision  tactique  offrant  des  previsions  exactes  et  coherentes  sur  les 
conditions  acoustiques  et  la  detectabilite  des  cibles.  Le  progiciel  Bellhop  est  au  cceur  du  systeme 
qui  offre  des  previsions  acoustiques  aux  programmes  clients. 

Perspectives 

II  est  prevu  de  poursuivre  les  ameliorations  du  programme  Bellhop  en  integrant  des  modeles  plus 
exacts  qui  seront  plus  efficaces  d’un  point  de  vue  operationnel. 


IV 
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1.  Introduction 


Bellhop  is  a  computer  program  created  by  Dr.  Michael  Porter  that  computes  acoustic  fields  in 
oceanic  environments  via  Gaussian  beam  tracing.  The  environment  consists  of  an  ocean  that  may 
have  range  variations  in  the  sound  speed  profile,  the  bottom  loss,  and  the  bathymetry.  Two 
programs  were  created  for  use  with  the  Environment  Modeling  Manager  (EMM)  of  the  Canadian 
Navy  at  the  Defence  R&D  Canada  -  Atlantic  (DRDC  Atlantic)  laboratory.  The  first  is  the  passive 
version  of  Bellhop  named  BellhopDRDC_ray_TL_v4.  The  outputs  from  the  passive  program 
include  transmission  loss  (coherent,  semi-coherent,  and  incoherent)  and  ray  traces.  The  second  is 
the  active  version  of  Bellhop  named  BellhopDRDC_active_v4.  The  outputs  from  the  active 
program  include  the  arrival  tables,  reverberation  time  series,  target  echo  time  series,  and  the 
signal  excess  versus  range. 

Changes  in  version  4: 

The  changes  between  this  fourth  version  and  the  previous  version  [1]  include: 

1.  The  addition  of  a  range  interpolation  scheme  for  range  dependent  sound  speed 
profiles. 

2.  The  addition  of  curvilinear  interpolation  for  bathymetry. 

3.  The  incorporation  of  changes  from  Dr.  Porter’s  Web  version  of  Bellhop  dated  May 
2008. 

The  Fortran  coding  approaches  are  similar  to  the  previous  versions.  They  consist  of  a  frontend 
program  that  reads  the  input  files  and  writes  the  output  files,  and  a  subroutine  named 
BellhopDRDC_*_v4  (where  *  represents  either  ‘active’  or  ‘ray_TL’).  This  structure  was  used  to 
enable  repeated  calls  to  the  subroutine  Bellhop  from  within  the  frontend,  for  looping  calculations 
over  source  depth,  frequency,  or  bearing,  for  example.  It  is  anticipated  that  the  user  may  rewrite 
or  replace  the  frontend  algorithm  to  suit  his  own  needs. 

Also  included  in  this  users  guide  is  a  simplified  program  to  compute  surface  and  bottom  loss  for 
separate  analysis.  Finally,  also  included  are  examples  of  simple  IDL  plot  routines  for:  TL  vs 
range,  (full  field  or  single  depth),  ray  tracing,  an  SSP  map,  reverberation,  signal  excess,  surface 
loss,  surface  scattering  strength,  bottom  loss,  and  bottom  scattering  strength. 
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2.  BellhopDRDC ray TL v4 


The  BellhopDRDC_ray_TL_v4  model  is  intended  for  passive  predictions  of  ray  paths  and 
transmission  loss.  This  model  consists  of  five  Fortran  source  files  and  their  subroutines: 

1 .  datamod_ray_TL_v4.f90  -  module  with  data  array  declarations 

2.  refcomod_ray_TL_v4.f90  -  module  with  reflection  coefficient  array  declarations  and  some 
loss  models 

CALCbotRC  -  compute  bottom  reflection  coefficients  using  MGS 
NAVOCEANO  routine 

CALCtopRC  -  compute  surface  reflection  coefficient  using  either  Modified 
Eckart  or  Beckmann-Spezzichino 
BOTT_NEW  -  MGS  bottom  loss  function 

SURF_NEW  -  Surface  loss  function  using  Bechmann-Spezzichino 
LFSOPN  -  Low  Frequency  Open  Ocean  Surface  loss  using  Modified  Eckart 

3.  sspmod_ray_TL_v4.f90  -  module  with  sound  speed  array  declarations 

4.  frontend_ray_TL_v4.f90 

frontend_ray_TL_v4  -  main  program  and  outputs  to  Bellhop.log  and  *TL.txt. 
clean_up  -  deallocate  ray  structure  arrays 
Ray  write  -  write  out  to  rays.txt  the  ray  trace  path  information 
READIN_v4  -  reads  file  runinput_v4.inp  and  allocates  and  initializes  arrays  for 
range  and  receiver  depth 

READBTY_v4  -  reads  file  bathy.inp  and  allocates  arrays  for  bathymetry 
READSVP_v4  -  reads  file  speed,  inp  and  allocates  arrays  for  sound  speed 
READBOTLOSS_v4  -  reads  file  bottomloss.inp  and  allocates  arrays  for  which 
ever  bottom  type  was  specified 

READBPATTERNS_v4  -  reads  file  becimpattem.inp  for  sensor  beam  pattern, 
allocates  arrays  and  converts  loss  to  pressure  coefficient 
TMP_SPP  -  function  to  convert  temperature  to  sound  speed  using  Leroy’ s 
equation 

dumpsspmap  -  write  out  to  SSPmap.txt  the  SSP  sampled  over  range  and  depth 
for  contour  plot 

5.  bellhopDRDC_ray_TL_v4.f90 

BellhopDRDC_ray_TL_v4a  -  beginning  of  bellhop  algorithm-  initializes  arrays, 
calls  ray  trace  and  calls  Grab  style  TL  computation 
Trace  -  traces  a  ray  for  each  launch  angle 
Step  -  takes  a  single  step  along  the  ray  path 
Reducestep  -  computes  step  size  to  land  on  key  points 

Reflect  -  changes  ray  direction  and  computes  amplitude  and  phase  at  reflection. 

The  two-layer  geoacoustic  bottom  loss  is  embedded  in  this  subroutine. 
REFCO  -  interpolates  for  reflection  coefficients  from  table  if  needed 
INFLUGRB  -  Gaussian  beam  contribution  to  complex  pressure  for  TL 
Quad  -  chooses  method  of  interpolation  of  sound  speed  with  depth  and  range 
Linear  -  preferred  method  of  SSP  range  interpolation 
Smoother  -  Savitsky-Golay  smoothing  filter  for  coherent  TL 
Thorpe  -  Thorpe  attenuation 

CRCI  -  converts  real  wave  speed  and  attenuation  to  a  single  complex  wave  speed 
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ERROUT  -  outputs  error  messages 

The  three  module  files,  datamod_ray_TL_v4.f90,  refcomod_ray_TL_v4.f90,  and 
sspmod_ray_TL_v4.f90  contain  the  data  arrays  and  declarations,  and  must  be  compiled  first.  The 
executable  is  named  BellhopDRDC_ray_TL_v4.exe 

To  run  the  program,  place  the  executable  BellhopDRDC_ray_TL_v4.exe  in  your  working 
directory  or  on  your  path.  Place  the  five  input  files  listed  below  in  your  working  directory.  Then 
click  on  the  .exe  icon  or  use  the  windows  start/run  command.  If  programming  in  IDL,  the  spawn 
command  can  be  used  to  run  the  executable.  For  example,  the  command  to  run  this  in  IDL  is: 
spawn,  'BellhopDRDC_ray_TL_v4.exe',  result,  /noshell. 

2.1  Input  files 

There  are  five  input  files:  runinput_v4.inp,  speed.inp,  bottomloss.inp,  bathy.inp  and 
beampattern.inp.  The  formats  are  free  field,  so  the  values  on  each  row  do  not  occupy  specific 
column  positions,  but  only  need  be  separated  by  a  space. 

2.1.1  Runinput_v4.inp 

This  file  contains  scenario  and  runtime  choices,  as  defined  in  Table  1.  In  this  table,  the  following 
alphabetic  choices  are  defined: 

‘X1X2X3X4’  =  the  run  choice  options,  consisting  of  4  letters: 

X 1  =  type  of  output 

C  =  Coherent  transmission  loss  in  output  file  CTL.  txt 
S  =  Semi-coherent  transmission  loss  in  output  file  STL.  txt 
I  =  Incoherent  transmission  loss  in  output  file  ITL.  txt 
R=  Ray  trace  path  information  in  output  file  rays.txt 
X2  =  SSP  range  interpolation  method 
N  =  none,  uses  abrupt  change 
L  =  Linear,  the  preferred  method  and  default 
X3  =  Bathymetry  range  interpolation  method 

L  =  piecewise  linear,  the  preferred  method  and  default 
C  =  curvilinear  interpolation 
X4  =  Flag  for  using  Thorpe  volume  attenuation 

T  =  use  Thorpe  attenuation,  the  preferred  choice  and  default 
N  =  use  no  attenuation 

‘S’  =  surface  loss  model  choice 

‘B’  =  Beckmann  Spezzichino  surface  loss 

‘E’  =  Modified  Eckart  low  frequency  open  ocean  surface  loss,  default  model 
(Note  that  the  bottom  loss  model  is  chosen  in  the  bottomloss.inp  file) 
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Table  1.  runinput_v4.inp  file  structure 


Line  #,  entry 

Notes 

1. 

title 

up  to  70  characters  enclosed  in  single  quotes 

2. 

frequency 

Hz 

3. 

source  depth 

Meters 

4. 

number  of  receiver  depths 

# 

5. 

top  and  bottom  of  receiver  depth  array 

Meters-  note:  needs  slash  at  end  value  to  denote  an  array-  a  single 
value  can  also  be  used 

6. 

Range  step  for  output;  longest  range 

Meters;  Kilometres 

7. 

wind  speed;  surface  loss  model  choice 

‘S' 

Knots;  ‘S’  ={B,E} 

8. 

Run  Choice  options  ‘X1X2X3X4’ 

‘{C,S,I,R}  {L,N}  {L,C}  {T,N}’  choosing  one  letter  from  each  group  to 
comprise  a  4  letter  sequence  in  single  quotes 

9. 

Internal  step  size;  number  of  rays;  start 
angle;  stop  angle;  kill-after-bounce 
number 

Default  value  =  -1  Internal  step  size  in  m;  angles  in  degrees;  negative 
angles  first.  Default  is  -15  to  15  deg,  and  100  bounces 

For  ray  tracing,  the  number  of  rays  and  start  and  stop  angle  should  be 
selected  by  the  user.  For  transmission  loss,  these  should  be  defaulted 
to  -1. 

10 

Range  smoothing  flag,  dumpSSP  flag 

Meters,  Default  =  -1 ,  no  smoothing  Smoothing  only  affects  the  ‘C‘ 
coherent  TL 

Dumpssp  flag  =  1  to  write  out  SSP  in  range  and  depth  into  the  file 
SSPMap.txt 
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2.1.2  Speed,  inp 

This  file  contains  sound  speed  profiles  in  depth  and  range. 


Table  2.  speed. inp  file  structure 


Line  #,  entry 

Notes 

1 

Number  of  range  dependent  profiles 

2 

Range  to  profile;  number  of  points  in  that  specific  profile,  n 

km 

3  to  n+3 

Depth;  speed  or  temperature 

M;  m/sec  or  QC 

Repeat  from  2  for  each  profile 

Note:  there  should  always  be  a  point  at  the  surface  and  at  or  below  the  deepest  bathymetry  point. 
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(» continued  from  previous  page) 
250.0  1490.6 


Figure  2.  Sample  speed. inp  file. 

The  speed  may  be  plotted  on  a  single  graph,  spaced  10  m/sec  apart  using  the  plot  routine 
read_plot_speed.pro.  For  this  sample  SSP  file,  the  result  is  shown  in  Figure  3. 


Figure  3.  SSP’s  plotted  from  speed,  inp. 


2.1.3  Bottomloss.inp 

This  file  contains  the  range  dependent  bottom  loss  descriptions. 


Table  3.  bottomloss.inp  file  structure 


Line  #,  entry 

Notes 

1 

Bottom  treatment  option;  attenuation  units  ‘XY’ 

‘XY’:  ‘X’=  {M,A,T}  ‘Y’=  {F,M,W,N} 

2 

number  of  range  dependent  bottom  sets,  n 

3  to  n+3 

If  ‘X’=‘M’:  range;  province  number 

If  ’X’='A’:  range;  cl ;  rhol ;  attenl ;  hi ;  c2;  rho2;  atten2 

lf‘X’=’T’:  range;  #  of  table  rows 

Angle;  reflection  coefficient;  phase 

Km;  MGS  province  number 

Km;  m/sec;  g/cc;  units  of  ‘Y’;  m;  m/sec; 
g/cc;  units  of  ‘Y’ 

Km;  number  of  rows 

Degrees;  decimal  fraction;  degrees 
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In  this  table,  the  following  are  defined: 

‘X’  =  the  bottom  treatment  option 

'M'  =  MGS  or  HFBL  provinces 

'A'  =  Two  Geoacoustic  fluid  layers  (no  shear) 

'T'  =  Read  in  table  of  pressure  reflection  coefficients  and  phases  as  a  function  of 
grazing  angle 

‘Y’  =  The  attenuation  units  which  are  used  in  the  geoacoustic  layers  only,  choices  are 
F'  =  dB/(m  kHz) 

'M'  =  dB/m 

'W'  =  dB/wavelength 

'N'  -  nepers/m 


‘M’ 

!  Bottom  option  for  MGS 

3 

!  number  of  range  dependent  bottom  provinces 

0.  4 

!  range  (km),  province  number 

10.  2 

50.  8 

'AF' 

!  A=geoacoustic,  F=  dB/m  kHz 

6 

Inumber  of  bottom  regions 

0.0  1453.  1.41 

0.038  10.  1557.  1.73  0.156 

7.4  1807.  2.175  0.131  30.  3500.  2.60  0.020 

14.8  1630.  2.00 

0.157  4.  1880.  2.175  0.085 

16.7  1772.  2.11 

0.142  10.  1880.  2.175  0.085 

44.4  1630.  2.00 

0.157  4.  1880.  2.175  0.085 

48.1  1807.  2.175  0.131  30.  3500.  2.60  0.020 

! range  cl  rhol 

attenl  depth  c2  rho2  atten2 

T 

!  Bottom  option  for  table  of  reflection  coefficients  vs  angle 

1 

!  number  of  sets  of  tables 

0.  5 

!  range  (km),  number  of  entries  in  table 

0.  1.0  0.0 

!  angle  (deg),  reflection  coef  fraction,  phase  (deg) 

10.0.8  0.0 

30.0.7  0.0 

50.0.5  0.0 

90.0.5  0.0 

Figure  4.  Three  samples  of  bottomloss.inp  files. 
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2.1.4  Bathy.inp 

This  file  contains  the  bathymetry. 


Table  4.  bathy.inp  file  structure 


Line  #,  entry 

Notes 

1 

Number  of  bathymetry  points,  n 

2  to  n+2 

Range;  depth 

Km;  m 

Note:  needs  a  point  at  zero  range 


2.1.5  Beampattern.inp 

This  file  contains  the  receiver  vertical  beam  pattern  in  dB. 


Table  5.  beampattern.inp  file  structure 


Line  #,  entry 

Notes 

1 

Number  of  vertical  angles,  n 

2  to  n+2 

Angle;  loss 

Deg;  dB 

3 

Inumber  of  angles 

-90.  0.0 

langle(deg);  loss(dB) 

0.  0.0 

90.  0.0 

Figure  6.  Sample  beampattern.inp  for  an  omni-directional  beam. 
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2.2  Output  files 


There  are  six  possible  output  files  from  BellhopDRDC_ray_TL_v4.  The  computed  data  is  written 
to  .txt  files  in  ASCII,  depending  on  the  runtime  choices  made  in  the  input  file  runinput_v4.inp. 

ITL.  txt  created  by  run  choice  ‘I’ 

STL. txt  created  by  run  choice  ‘S’ 

CTL.txt  created  by  run  choice  ‘C’ 
rays.txt  created  by  run  choice  ‘R’ 

SSPmap.txt  created  by  dumpsspflag  =  1 
bellhop.log 

2.2.1  CTL.txt,  ITL.txt  or  STL.txt 

This  file  contains  the  transmission  loss  (either  coherent,  semi-coherent  or  incoherent,  depending 
on  the  choice  made  in  runinput_v4.inp).  At  the  top,  it  lists  the  run  title,  frequency  and  source 
depth.  The  next  line  contains  the  number  of  ranges  and  number  of  receiver  depths.  Following  this 
are  listed  the  range  array  in  km,  then  the  receiver  depth  array  in  m,  then  transmission  loss  in  dB 
by  range  and  receiver  depth.  An  example  listing  is  shown  in  Figure  7. 


BELLHOP- 

Emerald  basin  toward  Sambro  Bank 

1 200Hz  2 l.O.m  source  depth 

200 

200 

0.2500000 

0.5000000 

0.7500000 

1.000000 

1.250000 

1.500000 

1.750000 

2.000000 

2.250000 

2.500000 

2.750000 

3.000000 

3.250000 

3.500000 

3.750000 

4.000000 

4.250000 

4.500000 

4.750000 

5.000000 

49.00000 

49.25000 

49.50000 

49.75000 

50.00000 

0.000000 

1.256281 

2.512563 

3.768844 

5.025126 

6.281407 

7.537688 

8.793970 

10.05025 

11.30653 

244.9749 

246.2312 

247.4874 

248.7437 

250.0000 

200.0000 

57.63440 

63.03972 

68.00800 

71.89261 

74.56378 

78.60585 

70.51897 

71.89424 

73.83934 

Figure  7.  Portion  of  an  ITL.txt  output. 


In  the  output  sometimes  the  first  several  transmission  loss  values  are  200dB,  as  shown  above. 
This  default  loss  occurs  if  the  first  depth  point  was  high  above  the  source  and  the  trace  angles 
were  defaulted  (runinput_v4.inp,  line  9,  start  and  stop  angle)  to  be  ±15°,  therefore  this  point  might 
not  have  been  ensonified  in  a  downward  refracting  profile.  The  same  default  loss  can  occur  at  a 
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deep  depth  point  below  the  source.  To  provide  very  short  range  loss  values  it  is  necessary  to  open 
up  the  angle  fan  to  ±25°  or  more,  at  the  cost  of  some  runtime. 


Bellhop  TL  vs  range 


Figure  8.  Example  transmission  loss  plot  for  60m  receiver.  Black  is  coherent,  CTL.txt.  Red  is 
semi-coherent,  STL.txt  and  blue  is  incoherent,  ITl.txt,  using  200  range  points. 


Figure  9.  Left:  example  of  full  field  plot  ofITL.txt  ( Incoherent  calculation)  which  was  computed 
using  200  receiver  depths  from  0  to  250m.  The  bathymetry  is  plotted  as  a  line  along 
the  bottom.  Right:  full  field  plot  ofCTL.txt  ( coherent  calculation).  Both  used  the 
linear  SSP  range  interpolation  and  linear  bathymetry  interpolation  with  a  70m  source 
depth. 
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Figure  9  displays  a  good  example  of  a  potential  pitfall  in  using  this  range  dependent  model.  Note 
that  in  the  figure,  there  are  places  where  the  field  extends  below  the  bathymetry,  since  the 
receiver  array  was  defined  to  250m  to  cover  the  deeper  part  of  the  water,  but  the  bathymetry  then 
rises  to  1 10m.  The  portion  of  the  field  below  the  bathymetry  is  not  a  true  representation  of  the 
acoustic  field  there.  The  loss  generated  by  Bellhop  on  reflection  from  the  bottom  into  the  water 
column  is  correct,  however  the  field  shown  within  the  bottom  does  not  have  the  right  level.  It  is 
an  artifact  of  the  Gaussian  beam  representation  in  Bellhop.  It  should  be  ignored  or  blanked  out  in 
the  figure,  and  in  all  other  Bellhop  applications,  care  should  be  exercised  that  the  user  is  only 
working  with  transmission  loss  values  from  those  receivers  positioned  above  the  bathymetry. 
When  receivers  are  defined  that  extend  below  the  bathymetry  at  some  point,  a  warning  is 
generated  and  written  to  bellhop.log. 

2.2.2  Rays.txt 

The  output  file  named  rays.txt  contains  ray  tracing  information.  Its  structure  is  to  echo  some  of 
input  choices  in  the  first  few  lines.  The  number  of  rays  being  traced  is  listed  (in  the  case  shown 
below  it  is  6).  Then  in  a  loop  over  the  number  of  rays,  each  ray  is  described  by  the  launch  angle  (- 
10.0)  and  number  of  steps  or  points  in  the  trace  (4521).  Finally,  the  [r,z]  coordinates,  ray  angle, 
delay  time,  and  number  of  surface  and  bottom  bounces  of  each  ray  are  listed  for  each  step.  Both  r 
and  z  are  given  in  m,  angle  is  in  degrees,  and  time  is  given  in  seconds. 

Figure  10  shows  a  portion  of  the  rays.txt  listing  for  the  70m  source,  and  it  demonstrates  an 
anomaly  that  always  occurs  in  Bellhop  ray  traces.  That  is  that  there  are  often  a  number  of 
repeated  points  (  see  for  example  the  line  at  184.1597m)  that  result  as  Bellhop  tries  to  place  a  ray 
exactly  on  a  sound  speed  depth  or  a  defined  bathymetry  range.  The  new  subroutine  in  Bellhop 
called  reducestep.f90  is  responsible.  It  does  not  affect  the  result  but  it  does  enlarge  the  file  sizes. 

The  rays.txt  output  can  be  plotted  with  the  bathymetry,  as  shown  in  Figure  1 1 .  The  case  shown 
was  computed  using  the  default  20  rays  from  -10°  to  +10°,  with  a  70m  source  depth  so  that  the 
figure  would  correspond  directly  to  the  full  field  transmission  loss  plot  in  Figure  9. 


BELLHOP- 

Emerald  basin  toward  Sambro  Bank 

1200.0Hz 

70.0m  source  depth 

Kill  Trace  after  50  bottom  bounces 

20 

-10.00000 

4735 

0.0000000E+00  70.00000 

-9.999999 

0.0000000E+00 

0  0 

1.47721 17E-07  70.00000 

-9.999999 

1.0218958E-10 

0  0 

14.77122 

67.39022 

-10.03921 

1.0219570E-02  0 

0 

29.54066 

64.77034 

-10.07839 

2.0440368E-02  0 

0 

184.1567 

40.00029 

-5.499722 

0.1269062  0 

0 

184.1597 

40.00000 

-5.499527 

0.1269083  0 

0 

184.1597 

40.00000 

-5.499527 

0.1269083  0 

0 

199.1020 

38.67991 

-4.599111 

0.1370111  0 

0 

.  .  . 

Figure  10.  Portion  of  a  rays.txt  output. 
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Bellhop  Ray  Plot 


Figure  11.  Plot  of  rays.txt  for  a  70m  source  showing  the  reflections  from  the  uneven  bathymetry. 

2.2.3  Bellhop.log 

This  file  contains  a  log  of  the  runtime  statements  generated  in  any  run.  Some  inputs  are  echoed, 
and  any  warnings  or  errors  are  listed  here  as  generated  by  the  Bellhop  code. 

BELLHOP-  Emerald  basin  toward  Sambro  Bank 

Frequency=  1200.000 
Source  depth=  70.00000 

range  step(m)  =  250.0000  Maximum  range(km)=  50.00000 
Wind  speed  (kts)=  10.00000 
Beckman-Spezzichino  surface  loss 
Runchoice=  Ray  trace 

Thorpe  volume  attenuation  used  for  frequency  dependent  water  column  absorption 
No  range  smoothing 

range, depth  computed  SSP  matrix  written  to  SSPmap.txt 
Number  of  receiver  depths=  200 

Top  and  bottom  of  Receiver  depths=  0.0000000E+00  250.0000 
Bathymetry  interpolation  is  linear  piecewize 
Number  of  bathymetry  points=  12 

Range(km)  Depth(m) 

1  0.0000000E+00  238.0000 

2  5.500000  256.0000 

3  18.50000  238.0000 

f>>  Continued  on  next  page) 
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(>>  Continued  from  previous  page) 

***  WARNING  *** 

Generated  by  program  or  subroutine:  Bathy.inp 
Receiver  deeper  than  bathymetry 

Number  of  sound  speed  profiles=  3 
Linear  range  interpolation  used  on  SSP 
Range(km)=  0.0000000E+00 

1  0.0000000E+00  1498.000 

2  30.00000  1499.200 

3  35.00000  1491.700 

Bottom  option=  Acoustic  parameters 
Atteunation  unit  choice=  dB/(m  kHz) 

Number  of  range  dependent  bottom  properties=  6 

Range(km)=  0.0000000E+00 
c2,rho2,a2,h2,c3,rho3,a3 

1453.000  1.410000  3.7999999E-02  10.00000  1557.000 

1.730000  0.1560000 

Range(km)=  7.400000 
c2,rho2,a2,h2,c3,rho3,a3 

1807.000  2.175000  0.1310000  30.00000  3500.000 

2.600000  2.0000000E-02 

Sensor  Beampattem 
angle(deg),  bpat(dB) 

-90.00000  0.0000000E+00 

0.0000000E+00  0.0000000E+00 
90.00000  0.0000000E+00 

Successful  input  read 

BELLHOP-  Emerald  basin  toward  Sambro  Bank 

Number  of  rays  =  20  from  -10.00000  deg  to  10.00000  deg 

Kill-after-bounce  50 
Minimum  Step  size(m)  =  15.0000000000000 

CPU  Time  =  1.07  seconds 


Figure  12.  Sample  portions  of  a  bellhop.log. 


2.2.4  Sspmap.txt 

This  file  contains  a  200x200  sample  map  of  the  sound  speed  profile  with  range  using  the 
interpolation  scheme  selected  in  runinput_v4.inp,  line  8.  The  first  line  of  sspmap.txt  lists  the 
number  of  range  and  depth  points  and  a  letter  indicating  the  type  of  range  interpolation,  L=linear 
and  N=none.  Then  the  range  points  are  listed,  followed  by  the  depth  points,  followed  by  the 
sound  speed  in  range  and  depth. 


DRDC  Atlantic  CR  201 0-1 34 


13 


200  200  L 

0.0000000E+00  0.2500000  0.5000000  0.7500000  1.000000 

1.250000  1.500000  1.750000  2.000000  2.250000 

2.500000  2.750000  3.000000  3.250000  3.500000 


Figure  13.  Sample  portion  of  the  sspmap.txt  file. 


Figure  14.  Example  of  sspmap.txt,  a  contour  plot  of  the  SSP  with  range  and  depth. 


2.3  Plot  routines 

Several  IDL  plot  routines  have  been  prepared  to  provide  a  simple  graphic  representation  of  the 
output  products  from  BellhopDRDC.  These  should  be  freely  altered  to  suit  the  users’  data  and 
output  requirements. 

Read  tl  plot  loss.pro:  Routine  to  read  each  of  the  xTL.txt  output  files,  as  in  the  example  in 
Figure  8.  The  user  will  be  asked  to  enter  the  receiver  depth.  If  it  does  not  exactly  match  one  of 
the  computed  depths,  the  plot  routine  will  choose  the  next  closest  depth.  Presently  the  plot 
routine  is  set  to  open  each  of  the  three  xTL.txt  files  and  over  plot  them  all  in  color. 
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Read  tl  plot  field.pro:  Routine  to  read  the  xTL.txt  output  and  the  bathy.inp  file  as  in  the 
examples  in  Figure  9.  Presently  the  plot  routine  is  set  to  open  each  of  the  three  xTL.txt  files  and 
plot  each  in  a  separate  window. 

Read  rays  plot  trace. pro:  Routine  to  read  the  output  in  rays.txt  and  the  bathymetry  in  bathy.inp 
and  produce  a  ray  trace  figure  as  shown  in  Figure  1 1 . 

Read  plot  speed.pro:  Routine  to  read  the  input  file  of  SSP,  speed.inp,  and  overplot  all  the 
profiles  spacing  them  lOm/s  apart  as  shown  in  Figure  3. 

Read  sspmap.pro:  Routine  to  read  the  sspmap.txt  that  was  created  if  the  dumpsspmap  option  was 
selected  in  runinput_v4.inp.  The  plot  shows  a  contour  map  of  the  SSP  in  range  and  depth  as 
shown  in  Figure  14  with  the  bathy.inp  file  overplotted. 
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3.  BellhopDRDC active v4 


The  BellhopDRDC_active_v4  model  is  intended  for  active  predictions  of  bistatic  target  echo  time 
series,  bistatic  reverberation  and  active  signal  excess  using  SALT  (Sound  Angle,  Level  and  Time) 
tables  produced  by  the  incoherent  output  from  Bellhop.  This  model  consists  of  fourteen  Fortran 
source  files  and  their  subroutines: 

1 .  datamod_active_v4.f90  -  module  of  data  array  declarations  and  size  limitations 

2.  refcomod_active_v4.f90  -  module  of  reflection  coefficient  array  declarations 

CALCbotRC  -  computes  bottom  reflection  coefficients  using  MGS 
NAVOCEANO  routine 

CALCtopRC  -  computes  surface  reflection  coefficient  using  either  Modified 
Eckart  or  Beckmann-Spezzichino 
BOTT_NEW  -  MGS  bottom  loss  function 

SURF_NEW  -  Surface  loss  function  using  Bechmann-Spezzichino 
LFSOPN  -  Low  Frequency  Open  Ocean  Surface  loss  using  Modified  Eckart 

3.  saltmod_active.f90  -  module  with  SALT  table  array  allocation  declarations 

4.  SEmod_active.f90  -  module  with  SE  input  variable  allocation  declarations 

5.  SSPmod_active_v4.f90  -  module  with  sound  speed  array  allocation  declarations 

6.  frontend_active_v4.f90 

frontend_active_v4  -  main  program 

Setdefaults  -  assigns  default  inputs  for  active  applications 

7.  readinput_active.f90 

readinput_active  -  reads  input  files  for  SE,  speed,  bathy,  beam  patterns 
READBOTLOSS_S  -  reads  bottom  loss  and  allocates  arrays  for  whichever 
bottom  type  was  specified 
READreverb  -  read  user  input  reverberation  table 
CALCreverb  -  Rough  estimate  of  reverb  in  dB  using  401og(t)  fall-off 

8.  bellhopDRDC_active_v4.f90 

BellhopDRDC_active_v4  -  beginning  of  bellhop  algorithm  which:  initializes 

arrays,  calls  ray  trace  and  calls  TL  computation,  and  defines  extra  receiver 
points  on  surface  and  bottom  (conforming  to  bathymetry)  for 
reverberation 

Trace  -  traces  a  ray  for  each  launch  angle 

Step  -  takes  a  single  step  along  the  ray  path 

Reducestep  -  refines  the  step  length  to  land  on  points  of  interest 

Reflect  -  changes  ray  direction  and  computes  amplitude  and  phase  at  reflection. 

The  geoacoustic  bottom  loss  is  embedded  in  this  subroutine 
REFCO  -  interpolates  for  reflection  coefficients  from  table  if  needed 
INFLUGRB  -  computes  Gaussian  beam  contribution  to  complex  pressure.  The 
point  on  the  bottom  is  shifted  at  each  step  to  conform  to  the  bathymetry. 
Results  are  sent  to  AddArr 
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QUAD  -  finds  sound  speed  and  gradient  using  interpolation  style  None  or  Linear 
Linear  -  Bilinear  quadrilateral  interpolation  of  SSP 

TMP_SPP  -  function  to  convert  temperature  to  sound  speed  using  Leroy’s 
equation 

Smoother  -  Savitsky-Golay  smoothing  filter 

AddArr  -  creates  arrival  SALT  table  for  surface,  bottom,  and  target  depths  from 
all  sensors  and  transmitter  along  all  bearings 
Thorpe  -  Thorpe  attenuation 

CRCI  -  converts  real  wave  speed  and  attenuation  to  a  single  complex  wave  speed 
ERROUT  -  outputs  error  messages 

9.  envstore_v4.f90  -  moves  range  dependent  environments  from  input  storage  arrays  into 
Bellhop  runtime  arrays  for  each  bearing  and  sensor.  Computes  internal  trace  step  size, 
deltas,  based  on  the  minimum  depth  of  the  bathymetry  on  that  bearing 

10.  reverb.f90  -  computes  bistatic  reverberation  from  surface  and  bottom  using  SALT  tables 
for  each  sensor  and  target  bearing.  Formulas  for  various  surface  and  bottom  scattering 
strengths  are  embedded.  Output  is  reverberation  time  series  without  source  level  for 
each  sensor 

1 1 .  scatstrength_v4  -  contains  all  scattering  strength  models 

OE  -  Ogden-Erskine  surface  scattering  strength 
CH  -  Chapman-Harris  surface  scattering  strength 
EC  -  Ellis-Crowe  bottom  scattering  strength 
LB  -  Lambert’s  rule  bottom  scattering  strength 
OM  -  Omni  bottom  scattering  strength 

12.  salt_v4.f90  -  stores  SALT  tables  for  each  sensor  and  bearing 

13.  SE_active.f90  -  computes  signal  excess  from  reverb,  target  echo  and  noise  for  each 
sensor  and  bearing.  Source  level  and  target  strength  are  applied.  The  result  is  saved  as  a 
function  of  range,  target  depth,  target  bearing  and  sensor. 

14.  targetecho.f90  -  computes  bistatic  signal  intensity  as  a  function  of  time,  target  range  and 
depth  along  target  bearing.  Output  is  signal  time  series  without  source  level  or  target 
strength. 

15.  writeoutput_active.f90 

WriteArrival  -  writes  SALT  arrival  tables  for  each  sensor.  Note  this  output  file 
is  only  a  portion  of  the  SALT  tables  on  the  target  bearing. 

Writerevb  -  writes  reverberation  time  series  for  each  target  bearing  and  sensor 
WriteSE  -  writes  SE  for  target  bearing,  target  depth,  range  and  sensor 
Writesignal  -  writes  target  echo  time  series  for  each  target  bearing,  target  depth, 
range  and  sensor 

WriteTL  -  writes  TL  from  transmitter  and  from  sensor  to  target  vs  range  for 
target  depth 

The  five  module  files,  datamod_active_v4.f90  and  refcomod_active_v4.f90,  saltmod_active.f90, 
SEmod_active.f90,  and  SSPmod_active_v4.f90  contain  the  data  arrays  and  declarations,  and  must 
be  compiled  first.  The  executable  is  named  BellhopDRDC_active_v4.exe 
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To  run  the  program,  place  the  executable  BellhopDRDC_active_v4.exe  in  your  working  directory 
or  on  your  path.  Place  the  five  input  files  listed  below  in  your  working  directory.  Then  click  on 
the  .exe  icon  or  use  the  windows  start/run  command.  If  programming  in  IDL,  the  spawn 
command  can  be  used  to  run  the  executable.  For  example,  the  command  to  run  this  in  IDL  is: 
spawn,  'BellhopDRDC_active_v4.exe',  result,  /noshell . 

3.1  Input  Files 

There  are  five  input  files:  active_general.inp,  radial_ssp.inp,  radial_bottomloss.inp, 
radial_bathy.inp  and  beampat_active.inp.  The  formats  are  free  field,  so  the  values  on  each  row  do 
not  occupy  specific  column  positions,  but  only  need  be  separated  by  a  space. 

For  active  use,  the  following  are  defaulted  in  the  file  frontend_active_v4,  subroutine  setdefaults: 
runchoice  =  T;  computes  incoherent  pressure 
Thorpe  =  ‘T’;  uses  Thoipe  attenuation 

numbotkill  =  100;  only  allow  up  to  100  surface  or  bottom  bounces 
angle  1,  angle2  =  ±25  deg;  range  of  up  and  down  angles  to  be  traced 
deltas0=  -1;  default  to  internally  calculate  the  ray  trace  range  step 
Nbeams0=  -1;  default  to  internally  calculate  the  number  of  rays  to  trace 

The  following  are  the  current  array  size  limitations  that  are  set  in  datamod_active_v4.f90: 

Nprofmax  =  25;  max  #  of  different  SSP’s  and/or  bottom  losses  along  any  single  bearing 

Nsspmax  =  200;  max  #  of  points  in  any  SSP 

NBathymax  =  500;  max  #  of  points  in  any  bathymetry  track 

Ntab  =181;  max  #  of  table  points  in  bottomloss  and  beampattern  table  input 

MxnArr  =  100;  max  #  of  arrivals  for  each  (depth.range)  in  SALT  tables 

Mxn  =  200000;  max  #  of  steps  in  each  ray  trace 


3.1.1  active_general.inp 

This  file  contains  the  basic  choices  for  the  scenario,  system  parameters,  scattering  strength 
models,  and  surface  loss  models. 

Table  6  lists  the  model  choices  and  Figure  15  shows  a  sampling  listing.  The  options  are: 

‘M’  =  surface  loss  model  choice 

‘B’  =  Beckmann  Spezzichino  surface  loss 

‘E’  =  Modified  Eckart  low  frequency  open  ocean  surface  loss  =  default  model 
(Note  that  the  bottom  loss  model  is  chosen  in  the  NUWbottomloss.inp  file) 

‘SM’  =  surface  scattering  strength  model  choice 

‘OE’  =  Ogden-Erskine  surface  scattering  strength-  a  combination  Chapman 
Harris  with  low  wind  speed  algorithms 

‘CH’  =  Chapman  Harris  surface  scattering  strength  =  default  model 
‘BM’  =  bottom  scattering  strength  model  choice 

‘EC’  =  Ellis  and  Crowe=  Lambert’s  rule  with  a  high  angle  facet  scattering  term= 
default  model 

‘LB’  =  Lambert’s  Rule  with  Mackenzie  Coefficient 
‘OM’  =  Omni-directional  Rule  with  Mackenzie  Coefficient 
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‘SB’  =  two  letters  for  the  interpolation  choices  for  ssp  and  bathymetry 
First  position  =  ssp  range  interpolation,  N=none,  L=linear 
Second  position  =  bathymetry  range  interpolation,  L=linear,  C=curvilinear 
Default  string  is  ‘LL’,  that  is,  both  interpolations  are  linear 


Table  6.  active _general.inp  file  structure 


Line  #,  entry 

Notes 

1 

Title 

80  characters  enclosed  in  single 
quotes 

2 

Number  of  receiving  sensors,  nsensor 

3 

Array  of  Noise  level  at  each  sensor 

Array  1  to  nsensor,  values  separated 
by  spaces,  units=dB 

4 

Source  level;  detection  threshold;  target  strength;  system  loss; 
pulse  length 

dB;  dB;  dB;  dB;  seconds 

5 

Array  of  Blast  arrival  time  at  each  sensor 

Array  1  to  nsensor,  values  separated 
by  spaces,  units=  seconds 

sign  will  position  the  sensor  to  the 
right  or  left  of  the  transmitter  as  you 
face  the  target  bearing.  Negative=left, 
positive=right 

6 

Frequency 

Hz 

7 

Array  of  Asset  depths  (sensors  and  transmitter) 

Array  1  to  nsensor  +1,  values 
separated  by  spaces,  units=m,  extra 
last  point  is  the  transmitter  depth 

8 

Maximum  range  to  target 

Km;  Note:  This  will  be  increased 
internally  to  include  the  max  distance 
between  sensor  and  transmitter  plus  a 
pulse  length. 

9 

Number  of  target  depths;  Target  depth  minimum;  Target  depth 
maximum 

Program  will  create  an  array  of  target 
depths,  units=  m,  Note:  must  have 
both  min  and  max  depth,  even  if  the 
number  of  target  depths=1 

10 

Wind  speed;  surface  loss  model  choice  M' 

kts;  ‘M'={B,E},  default=E 

11 

Surface  scattering  strength  model  choice  ‘SM’ 

‘SM’={OE,  CH},  default=CH 

12 

Bottom  scattering  strength  model  choice  ‘BM’;  Mackenzie 
coefficient;  Normal  incidence  bottom  loss;  facet  width  or  RMS  slope 

‘BM’={EC,  LB,  OM},  default=’EC’; 
dB,  used  with  all  choices; 
dB,  used  with  ‘EC’; 
degrees,  used  with  ‘EC’ 

13 

Range  interpolation  choices  for  SSP  and  Bathymetry  ‘SB’ 

ssp  =  { N , L} ;  Bathy  =  {L,C};  input  both 
choices  as  two  letter  string;  default  is 
‘LL’ 
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‘Test  run  with  1  sensor,  4  bearings’ 

1 

! Number  of  sensors 

70.8 

INoise  level  by  sensor,  dB 

230.0  0.0  0.0  0.0 

1.0  !SL,  DT,  TS,  Syslos,  in  dB,  pulse  length  in  sec 

-12.34 

!  blasttime  to  sensor,  right=+  left=- 

1200. 

!freq 

18.3  21. 

!  asset  depth  array,  including  transmitter 

50. 

!  maximum  range  to  target  locations  km 

2  21.0  80.0 

Inumber  of  target  depths,  dmin,  dmax 

10.  'E' 

Iwind  speed  in  knots,  surface  loss  model 

'OE' 

!  Surface  scattering  strength  model 

EC'  27.  5.  10. 

IBottom  scattering  strength  model  and  inputs  to  ‘EC’  model 

‘LL’ 

ISSP  interp  =  {N,L},  Bathy  interp  =  {L,C} 

Figure  15.  Example  of  active _general.inp  file. 


3.1.2  radial_ssp.inp 

This  file  contains  the  bearing  and  range  dependent  sound  speed  profiles  for  each  asset.  Currently 
the  dimensions  of  the  sound  speed  arrays  are  limited  to  at  most  25  profiles  along  each  bearing, 
each  profile  having  at  most  200  points.  At  least  one  of  the  profiles  along  a  bearing  track  must 
have  as  SSP  point  deeper  than  the  deepest  point  in  the  bathymetry  along  that  track.  Any  other 
profiles  on  that  track  will  be  linearly  interpolated  to  that  depth.  The  format  is  shown  in  Table  7 
and  a  sample  listing  is  illustrated  in  Figure  16. 


Table  7.  radial _ssp.inp  file  structure. 


Line  #,  entry 

Notes 

(no  input 
required) 

Loop  over  assets  (transmitter  last) 

Ensure  sensors  match  the  order  used 
in  radial_bathy.inp.  The  number  of 
recievers  nsensor  is  on  line  1  in 
active_general.inp.  The  number  of 
assets  is  nsensor+1  to  include  the 
transmitter 

0 

(no  input 
required) 

Loop  over  bearing  for  each  asset 

Note:  there  must  be  the  same  number 
of  bearings  for  all  assets.  The  number 
of  bearings  is  on  line  1  in 
radial_bathy.inp 

1 

Number  of  profiles  along  each  bearing 

Limited  to  25 

2 

Range  to  profile;  number  of  svp  points  n 

Km;  #  pts  limited  to  200 

3  to  3+n 

Depth;  speed  or  temperature 

m;  m/sec  or  degrees  C 

Repeat  from  line  1  for  next  bearing 

Repeat  from  line  0  for  next  asset 
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3 

!  first 

receiver,  first  bearing,  number  of  ssp  on  this  bearing 

0. 

18 

!  range  to  ssp  (km);  number  of  points  in  svp 

0.  1498.0 

30.0 

1499.2 

35.0 

1491.7 

40.0 

1483.8 

45.0 

1475.4 

50.0 

1466.5 

75.0 

1468.2 

80.0 

1470.0 

90.0 

1473.5 

95.0 

1475.3 

100.0 

1477.0 

125.0 

1479.6 

150.0 

1482.1 

175.0 

1484.6 

200.0 

1488.4 

225.0 

1489.5 

250.0 

1490.6 

300.0 

1490.5 

10.  5 

0. 

1499.0 

50. 

1468.0 

100. 

1470.0 

200. 

1485.0 

300. 

1491.0 

40.9 

0.  1498.0 

!depth(m)  speed(m/sec) 

35.0 

1491.7 

45.0 

1475.4 

75.0 

1468.2 

90.0 

1473.5 

100.0 

1477.0 

150.0 

1482.1 

200.0 

1488.4 

250.0 

1490.6 

1 

Inumber  of  ssp  on  this  bearing 

0. 

18 

190  deg  bearing,  18  points 

0.  1498.0 

30.0 

1499.2 

35.0 

1491.7 

Repeat  for  each  bearing  and  asset,  transmitter  last 

Figure  16.  Example  of  portion  of  radial _ssp.inp  file. 
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3.1.3  radial_bottomloss.inp 

This  file  contains  bottom  loss  information  for  each  asset  and  bearing.  The  type  of  bottom  selected 
(MGS,  geoacoustic,  or  table)  will  apply  to  all  assets  and  bearings.  Currently  the  dimensions  of  the 
loss  arrays  are  limited  to  25  different  regions  along  each  bearing,  and  for  table  entries,  the  number 
of  points  is  limited  to  91.  Table  8  shows  the  available  options  and  a  sample  listing  is  illustrated  in 
Figure  17.  The  options  are  : 

‘X’  =  the  bottom  treatment  option 

'M'  =  MGS  or  HFBL  provinces 

'A'  =  Geoacoustic  fluid  layers  (no  shear) 

'T'  =  Read  in  table  of  pressure  reflection  coefficients  and  phases  as  a  function  of 
grazing  angle 

‘Y’  =  the  attenuation  units  that  are  used  in  the  geoacoustic  layers  only,  choices  are: 

F'  =  dB/(m  kHz) 

'M'  =  dB/m 

'W'  =  dB/wavelength 

'N'  =  nepers/m 


Table  8.  radial Jbottomloss.inp  file  structure. 


Line  #,  entry 

Notes 

1 

Bottom  treatment  option;  attenuation  units 

‘XT:  ‘X’=  {M,A,T}  ‘Y’=  {F,M,W,N} 

(no  input 
required) 

Loop  over  assets  (transmitter  last) 

Ensure  sensors  match  the  order  used  in 
radial_bathy.inp.  The  number  of 
receivers  nsensor  is  on  line  1  in 
active_general.inp.  The  number  of 
assets  is  nsensor+1  to  include  the 
transmitter 

0 

(no  input 
required) 

Loop  over  bearing  for  each  asset 

Note:  there  must  be  the  same  number  of 
bearings  for  all  assets.  The  number  of 
bearings  is  on  line  1  in  radial_bathy.inp. 

2 

number  of  range  dependent  bottom  sets,  n 

Currently  limited  to  25 

3  to  n+3 

If  ‘X’=‘M’:  range;  province  number 

If  ’X’=’A’:  range;  cl ;  rhol ;  attenl ;  hi ;  c2;  rho2;  atten2 

If  ‘X’=T:  range;  #  of  table  rows;  then  loop  over  #  table 
rows  with  angle;  reflection  coefficient;  phase 

Km;  MGS  province  number 

Km;  m/sec;  g/cc;  units  of  ‘Y’;  m;  m/sec; 
g/cc;  units  of  ‘Y’ 

Km;  number  of  rows;  Degrees;  decimal 
fraction;  degrees 

Repeat  from  line  2  for  next  bearing 

Repeat  from  line  0  for  next  asset 
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!  attenuation  units  db/m  khz 

!0  deg  bearing  -  receiver  #1 


0.0  1453. 
26.8  1547. 
31.5  1630. 
40.7  1772. 
4 


1.41  0.038  10.  1557.  1.73  0.156 
1.72  0.158  20.  1880.  2.175  0.085 
2.00  0.157  4.  1880.  2.175  0.085 
2.11  0.142  10.  1880.  2.175  0.085 


!90  deg  receiver 


0.0  1453. 
24.0  1547. 
27.8  1807. 
46.3  1630. 
6 


1.41  0.038  10.  1557.  1.73  0.156 
1.72  0.158  20.  1880.  2.175  0.085 
2.175  0.131  30.  3500.  2.60  0.020 
2.00  0.157  4.  1880.  2.175  0.085 


!  1 80  deg  receiver 


0.0  1453.  1.41  0.038  10.  1557.  1.73  0.156 

14.8  1807.  2.175  0.131  30.  3500.  2.60  0.020 

19.4  1453.  1.41  0.038  10.  1557.  1.73  0.156 

33.3  1807.  2.175  0.131  30.  3500.  2.60  0.020 

44.4  1630.  2.00  0.157  4.  1880.  2.175  0.085 

50.0  1772.  2.11  0.142  10.  1880.  2.175  0.085 
6  !290  deg  receiver 

0.0  1453.  1.41  0.038  10.  1557.  1.73  0.156 


Figure  17.  Example  of  portion  of  radial  Jbottomloss.inp  showing  range  dependent  geoacoustic 
parameters  for  several  bearings. 


3.1.4  radial_bathy.inp 

This  file  contains  the  radials  and  bathymetry  for  all  radials  desired.  The  current  maximum  number 
of  bathymetry  points  on  any  radial  is  500.  A  flag  is  defined  to  mark  the  radial  bearing  which 
could  contain  the  target.  This  flag  will  be  used  to  trigger  the  computation  of  the  reverberation, 
target  echo  and  SE,  and  when  outputting  the  arrival  tables.  Table  9  shows  the  available  options 
and  a  sample  listing  is  illustrated  in  Figure  18. 
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Table  9.  radialjbathy.inp  file  structure. 


Line  #,  entry 

Notes 

1 

Number  of  radial  bearings 

Number  applies  to  all  sensors  and 
transmitter 

(no  input 
required) 

Loop  over  assets  (transmitter  last) 

The  number  of  receivers,  nsensor,  is 
on  line  1  in  active_general.inp.  The 
number  of  assets  is  nsensor+1  to 
include  the  transmitter 

0 

(no  input 
required) 

Loop  over  bearings  for  each  asset 

Note:  there  are  the  same  number  of 
bearings  for  each  asset 

2 

Radial  bearing;  number  of  bathymetry  points  n;  target  bearing  flag 

Degrees  measured  from  line  between 
sensor  and  transmitter;  #  of  points 
currently  limited  to  500  per  radial;  flag 
identifying  expected  target  bearing 

3  to  3+n 

Range;  depth  of  bathymetry 

Km;  m 

Repeat  from  line  2  for  next  bearing 

Repeat  from  line  0  for  next  asset 

4 

Inumber  of  radials-  same  for  all  assets 

0.  11 

0 

Hoop  over  assets,  radial  direction  phi(deg);#pts;  0/1  for  target  bearing 

0.0 

238 

1  range;  depth  of  bathy  for  #pts 

4.6 

256 

18.5 

238 

22.2 

219 

33.3 

183 

36.1 

165 

38.0 

146 

40.7 

128 

41.7 

110 

51.9 

91 

55.6 

91 

90.7 

0 

0.0 

238 

22.2 

219 

24.1 

183 

27.8 

165 

31.5 

146 

46.3 

128 

(>>  Continued  on  next  page) 
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3.1.5  beampat_active.inp 

This  file  contains  the  beampatterns  for  sensors  and  transmitter.  In  the  program,  these  beams  will 
be  assumed  to  be  pointing  along  the  target  bearings,  defined  by  the  target  bearing  flag  in 
radial_bathy.inp,  line  1.  Table  10  shows  the  available  options  and  a  sample  listing  is  illustrated  in 
Figure  19. 

There  is  presently  no  ability  to  specify  towed  array  beams. 


Table  10.  beampat_active.inp  file  structure. 


Line  #,  entry 

Notes 

(no  input 
required) 

Loop  over  assets  (transmitter  last) 

Ensure  sensors  match  the  order  used 
in  radial_bathy.inp.  The  number  of 
recievers  nsensor  is  on  line  1  in 
active_general.inp.  The  number  of 
assets  is  nsensor+1  to  include  the 
transmitter 

1 

Number  of  vertical  D/E  angles  in  pattern,  n 

(no  input 
required) 

Loop  over  number  of  angles  given  in  line  1 

2  to  2+n 

D/E  angle;  loss 

Degrees;  dB 

Repeat  from  line  1  for  next  asset 

DRDC  Atlantic  CR  201 0-1 34 
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3.2 

Output  files 

3.2.1 

Arrival.txt 

This  file  contains  incoherent  ray  arrival  structures  also  called  the  SALT  tables.  While  the  SALT 
tables  are  computed  for  all  bearings  and  assets,  this  output  is  only  triggered  on  the  target  bearing 
defined  by  the  flag  in  the  input  file  radial_bathy.inp.  For  a  configuration  of  BellhopDRDC  which 
outputs  all  the  SALT  tables,  this  trigger  should  be  set  to  1  on  all  bearings. 

As  illustrated  in  Figure  20,  this  file  begins  with  the  title  as  input  from  active_general.inp.  Next,  it 
lists  the  frequency  and  sensor  depth;  then  the  bearing  angle  and  sensor  number.  Next  it  lists  by 
column  the  target  depth(m),  range(km),  acoustic  intensity,  phase(rad),  delay  time(sec),  source 
angle(deg),  target  angle(deg),  number  of  reflections  from  the  surface,  and  the  number  of 
reflections  from  the  bottom.  A  header  with  abbreviations  of  these  outputs  is  given  for 
the  reader’s  convenience.  The  Fortran  output  format  for  these  numbers  is 


26 


DRDC  Atlantic  CR  201 0-1 34 


(f7.1,f7.2,2el2.4,f7.3,2f7.2,2i4).  This  listing  is  repeated  for  each  bearing  that  was  designated  a 
target  bearing  in  the  input  file  radial_bathy.inp  and  for  each  sensor  and  the  transmitter  (last 
listings).  The  write  statements  are  in  subroutine  Writearrival  in  the  file  writeoutput_active.f90. 
An  example  plot  of  the  arrival  angle  vs  range  is  plotted  in  Figure  2 1 . 

The  tables  include  all  the  target  depths  that  were  specified  in  active_general.inp  on  line  9.  They 
also  include  entries  for  the  surface  and  bottom  that  are  required  to  compute  reverberation.  The 
surface  entries  are  listed  first  at  0  depth,  then  come  the  target  depth  entries,  followed  by  the 
bottom  entries.  Note  that  the  depths  listed  for  the  bottom  entries  change  with  range  as  the  bottom 
contour  changes. 

The  amplitude  in  this  table  is  the  incoherent  acoustic  intensity.  The  transmission  loss  from  these 
entries  is  TL=10*log(sum  of  entries  at  the  same  range  and  depth).  It  is  possible  to  change  the 
default  run  choice  setting  to  produce  a  coherent  SALT  table,  however  in  that  case,  the  amplitude 
in  the  table  is  an  acoustic  pressure  and  the  phase  and  delay  time  must  be  used  to  produce  a 
coherent  intensity.  In  that  case,  the  computer  codes  in  the  other  active  products  of  reverberation, 
target  echo  time  series  and  signal  excess  that  are  programmed  to  use  intensity  inputs  would  also 
have  to  be  changed  to  a  coherent  calculation. 


‘Test  run  with  1  sensor,  4  bearings’ 

1 200Hz  1 8 .30m  source  depth 

290.  deg  Bearing  1  sensor  number 

Tdepth  Range  Intensity  Phase -rad  Time  Sangle  Rangle  Ntop  Nbot 

0.0  0.75  0.7292E-06  0.3142E+01  0.502  -0.89  1.96  1  0 

0.0  0.75  0.7174E-06  0.0000E+00  0.502  -0.20  -1.35  0  0 

0.0  0.75  0.1607E-07  0.0000E+00  0.502  0.10  -1.05  0  0 

0.0  0.75  0.7408E-09  0.0000E+00  0.502  0.30  -0.86  0  0 

0.0  0.75  0.7176E-09  0.0000E+00  0.498  1.58  10.12  0  0 

0.0  1.51  0.1 182E-09  0.5528E+01  1.067-16.60  16.99  2  1 
0.0  1.51  0.288 IE-08  0.5756E+01  1.067-16.40  16.81  2  1 


Figure  20.  Portion  of  output  file  arrival.txt  showing  some  of  the  surface  entries  using  Figure  15 
input. 
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Figure  21.  Example  of  arrival  angle  vs  range  plotted  using  arrival.txt  for  the  21m  target  depth. 


3.2.2  Reverb.txt 

This  file  contains  the  surface  and  bottom  reverberation  as  a  function  of  time,  with  the  source  level 
removed.  The  first  line  contains  the  title  from  active_general.inp.  The  next  line  gives  the  number 
of  time  points,  the  sensor  number  and  depth,  and  the  target  bearing  as  specified  in 
radial_bathy.inp.  Next  is  a  listing  of  the  time  array.  Then  follows  the  bottom  reverberation  time 
series  in  dB  (without  SL),  and  the  surface  reverberation  time  series  in  dB  (without  SL).  This  file 
structure  is  repeated  for  each  target  bearing  and  each  receiving  sensor.  The  write  statements  are  in 
subroutine  Writerevb  in  the  file  writeoutput_active.f90. 


The  -300.0000  value  shown  at  the  start  of  the  reverberation  section  of  Figure  22  is  a  default  value 


‘Test  run  with  1  sensor,  4  bearings’ 

189 

1  18.30000 

21.00000 

290.0000 

0.5000000 

1.000000 

1.500000 

2.000000 

2.500000 

3.000000 

3.500000 

4.000000 

4.500000 

5.000000 

5.500000 

6.000000 

6.500000 

7.000000 

7.500000 

-300.0000 

-300.0000 

-300.0000 

-300.0000 

-300.0000 

-156.5099 

-130.0203 

-128.0659 

-130.6896 

-131.9055 

-132.8520 

-136.2401 

-137.7021 

-140.5639 

-145.7615 

-147.2074 

-148.5177 

-150.4951 

-152.6780 

-154.4710 

Figure  22.  Selected  portions  of  Reverb.txt  output  file  using  Figure  15  loss  models. 
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designating  no  reception  at  that  time.  To  obtain  values  at  earlier  times,  the  defaulted  start  and 
stop  trace  angles  should  be  increased  from  ±25  deg  to  perhaps  ±45  deg. 

While  the  reverberation  is  an  important  computation  in  its  own  right,  this  file  is  primarily 
intended  to  be  combined  with  the  signal.txt  file  to  produce  the  signal  excess  that  may  be 
computed  at  a  later  date  or  in  another  language.  The  source  level  is  not  included  to  make  the 
signal  excess  computation  more  flexible. 

In  Figure  23,  the  top  plot  used  the  Beckmann  Spezzichino  ‘B’  surface  loss,  the  Chapman  Harris 
‘CH’  (surface)  and  Omni  ‘OM’  (bottom)  scattering  strengths  as  models  chosen  in 
active_general.inp.  By  way  of  contrast,  the  bottom  plot  used  Modified  Eckart  ‘E’  surface  loss, 
and  the  Ogden  Erskine  ‘OE’  (surface)  and  Ellis  Crowe  ‘EC’  (bottom)  scattering  strengths.  This 


1  sensor  at  depth=  18.3000 


0  5  10  15  20  25  30  35  40  45  50  55  60  65  70  75  SO  85  90  95 

iime(sec) 

Surfacefred-}  ana  BoitomCblaclO 


1  sensor  at  depth=  18.3000 


0  5  10  15  20  25  30  35  40  45  50  55  60  65  70  75  SO  85  90  95  100 

■t]me(sec) 

Surface  ('red')  and  Boito  black") 


Figure  23.  Plots  of  surface  and  bottom  reverberation  from  reverb.txt  (note:  source  level  is  not 
applied ).  The  top  plot  used  model  choices  {B,CH,OM}  while  the  bottom  plot  used 
{E,OE,EC . 
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latter  choice  of  models  produces  a  similar  contribution  from  the  knoll  bathymetry  of  the  bottom 
reverberation,  while  the  surface  reverberation  is  uniformly  higher. 


3.2.3  Signal.txt 

This  file  contains  the  target  echo  time  series  for  each  sensor’s  range,  target  depth  and  target 
bearing.  The  source  level  and  target  strength  are  not  included  in  this  output.  At  present  there  are 
no  graphics  that  make  use  of  this  file’s  output.  It  is  provided  to  be  an  input  along  with  the 
reverb.txt  file  for  signal  excess  calculations  that  may  wish  to  be  computed  at  a  later  date  or  in 
another  language.  The  source  level  and  target  strength  are  not  applied  to  this  file  to  make  the 
signal  excess  computation  more  flexible.  The  write  statements  are  in  subroutine  Writesignal  in 
the  file  writeoutput_active.f90. 

The  output  begins  with  the  run  title  from  active_general.inp.  Next,  looping  over  each  bearing  of 
interest  and  asset,  the  number  of  time  points,  the  number  of  ranges,  the  number  of  target  depths, 
the  sensor  number  and  target  bearing  are  listed,  followed  by  the  target  echo  time  series  in  dB. 


3.2.4  SE.txt 

This  file  contains  the  signal  excess  computed  using  the  Fortran  file  SE_active.f90  which  is 
included  in  this  program.  The  SE  computation  begins  by  working  in  intensity  units.  It  forms  the 
signal  term  by  multiplying  the  target  echo  time  series  by  the  source  level  and  target  strength.  It 
forms  the  interference  level  power  sum  by  adding  the  surface  and  bottom  reverberation  intensity 
time  series,  multiplying  by  the  source  level  and  adding  the  noise  intensity.  From  these  two  terms 
the  signal-to-interference  ratio  time  series  is  formed.  Taking  the  maximum  signal  to  interference 
ratio  from  the  time  series,  the  program  converts  this  to  dB  and  subtracts  the  dB  values  of 
detection  threshold  and  system  loss  to  form  the  signal  excess. 

The  SE  is  listed  in  the  file  SE.txt.  The  first  line  contains  the  title.  The  second  line  lists  the 
frequency  and  transmitter  depth.  The  next  line  contains  the  number  of  ranges,  the  number  of 
target  depths,  the  number  of  target  bearings  (those  that  were  identified  in  radial_bathy.inp),  and 
the  number  of  receiving  sensors.  The  range  array  is  listed  next.  Following  this,  the  output  loops 
over  sensor  number  and  writes  the  sensor  number  and  sensor  depth.  Inside  this  sensor  loop,  the 
output  loops  over  the  target  bearing  and  writes  its  value  in  degrees.  Lastly,  the  output  loops  over 
target  depth  and  writes  each  depth  in  meters,  followed  by  the  SE(dB)  vs  range  array.  The  write 
statements  are  in  subroutine  WriteSE  in  the  file  writeoutput_active.f90. 


‘Test  run  with  1  sensor,  4  bearings’ 


Frequency=  1200  transmitter  depth  (m)=  21.00000 
93  2  1  1 

0.7526418  1.505284  2.257925  3.010567  3.763209 

4.515851  5.268493  6.021134  6.773776  7.526418 


(>>  Continued  on  next  page) 
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(>>  Continued  from  previous  page) 

8.279059  9.031701  9.784343 

10.53699 

11.28963 

12.04227 

12.79491 

13.54755 

14.30019 

15.05284 

15.80548 

16.55812 

17.31076 

18.06340 

18.81605 

1  18.30000  sensor#  and  depth 

290.0000  target  bearing  (deg) 

2 1 .00000  target  depth  (m) 

-7.8558557E-02  -1.335800  2.14571 10E-02  0.7401745  2.898602 

8.519861 

6.849106 

4.837643 

9.012436 

11.73837 

9.486148 

10.76469 

12.86479 

7.269334 

9.710820 

13.33085 

11.35298 

7.134454 

9.388466 

10.77402 

10.35002 

10.18052 

8.118537 

9.292543 

10.42265 

80.00000 

-35.57817 

target  depth 
-37.65511 

(m) 

-37.85839 

-39.54664 

-36.73313 

-37.13431 

-38.05925 

-39.03345 

-40.48875 

-42.64893 

-42.27207 

-42.21347 

-43.35333 

-43.96018 

-45.02712 

-43.42149 

-45.56827 

-51.34261 

-48.63745 

-47.65379 

-49.21247 

-46.95230 

-51.52979 

-49.59116 

-50.81988 

Figure  24.  Example  of  portion  ofSE.txt  output. 


Figure  25.  Example  of  signal  excess  plot,  with  DT=TS-syslos-0.  Two  curves  are  for  two 
different  target  depths. 
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3.2.5  TL.txt 


This  file  contains  the  transmission  loss  as  a  function  of  range  for  the  transmitter  location  to  the 
target  depth,  and  for  the  sensor  location  to  the  target  depth  along  the  target  bearing.  The  beam 
patterns  of  the  transmitter  and  sensor  are  included  in  the  calculation.  Note  that  the  maximum 
range  is  larger  than  the  value  specified  in  the  active_general.inp  file.  Inside  the 
frontend_active_v4.f90  file,  the  maximum  range  is  automatically  increased  by  the  length  of  the 
baseline  between  sensor  and  transmitter  plus  a  pulse  length.  The  output  file  consists  of  the  title  of 
the  run,  followed  by  a  loop  over  sensor  number  and  target  bearing  in  which  the  next  two  lines 
contain  the  number  of  ranges,  the  number  of  target  depths,  then  the  sensor  number,  sensor  depth, 
transmitter  depth  and  the  target  bearing.  Following  this  is  the  range  array.  Then  in  a  loop  over 
target  depth,  the  depth  is  listed,  followed  by  the  transmitter  transmission  loss  array  vs  range,  and 
lastly  the  sensor  transmission  loss  array  vs  range.  The  write  statements  are  in  subroutine  WriteTL 
in  the  file  writeoutput_active.f90. 


‘Test  run  with  1  sensor,  4  bearings’ 


93 

1  18 
0.7526418 
4.515851 
8.279059 

68.49041 

21.00000 

55.33667 

63.31894 

67.26232 


2 

30000  21 

1.505284 
5.268493 
9.031701 

69.24305 

59.86486 

63.20134 

71.43450 


.00000  290 

2.257925 
6.021134 
9.784343 

69.99569 

59.65900 

67.48380 

67.33181 


.0000 

3.010567 

6.773776 

10.53699 


60.15039 

68.84501 

67.46899 


3.763209 

7.526418 

11.28963 


65.99412 

65.57792 

70.77728 


Figure  26.  Example  of  portion  of  TL.txt  output 
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transmitter  at  depth  =  21.0000 


Figure  27.  Example  ofTL  plot  from  the  TL.txtfile.  Left:  transmission  loss  from  transmitter  to 

target,  including  transmitter  beam  pattern  for  two  target  depths.  Right:  transmission 
loss  from  receiver  to  target,  including  receiver  beam  pattern  for  two  target  depths. 


3.2.6  Bellhop_active.log 

This  file  contains  a  log  of  the  runtime  statements  generated  in  any  run.  Some  inputs  are  echoed, 
and  any  warnings  or  errors  are  listed  here  as  generated  by  the  Bellhop  code. 
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3.3  Plot  Routines 


Some  IDL  plot  routines  have  been  prepared  to  provide  simple  graphic  representations  of  the 
active  program  outputs.  These  should  be  freely  altered  to  suit  the  users’  data  and  output 
requirements. 


Read_arrivals_plot_angleortimes.pro  -  reads  SALT  table  output  in  arrival.txt  and  plots  arrival 
angle  vs  range  as  shown  in  Figure  2 1 . 

Read_reverb_plot_re  verb. pro  -  reads  reverb.txt  and  plots  surface  and  bottom  reverberation  as 
shown  in  Figure  23. 

Read_SE_plot_SE.pro  -  reads  SE.txt  and  plots  signal  excess  vs  range  as  shown  in  Figure  25. 

Read_TL_plot_TL_active.pro  -  reads  TL.txt  from  the  active  program  output  and  plots  receiver  to 
target  transmission  loss  and  transmitter  to  target  transmission  loss  in  two  windows  as  a  function 
of  range  as  shown  in  Figure  27. 
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4.  Boundary  loss 


To  help  understand  the  transmission  loss  predictions  a  separate  program  called  boundary  loss  is 
provided  that  will  compute  the  bottom  and  surface  losses  and  the  surface  and  bottom  scattering 
strengths.  There  are  two  versions  of  this  program  for  active  and  passive. 

4.1  Boundaryloss_passive 

Boundaryloss_Passive.f90  -  main  program 
READIN_V4-  reads  runinput.inp 
READSVP_v4  -  reads  speed.inp 
READBOTLOSS_v4  -  reads  bottomloss.inp 

BOTT_NEW  -  function  to  compute  MGS  bottom  loss  from  province  numbers 
SURF_NEW  -  function  to  compute  surface  reflection  coefficients  using 
Bechmann-Spezzichino  formulas 

LFSOPN  -  computes  surface  loss  per  bounce  for  open  ocean  using  Modified 
Echart 

TMP_SSP  -  function  to  convert  temperature  to  speed  using  Leroy’ s  equation 
twolayerRefl  -  computes  reflection  coefficient  from  two  fluid  layers  of  sediment 
CRCIS  -  function  to  convert  real  wave  speed  and  attenuation  to  a  single  complex 
wave  speed 

ERROUT  -  outputs  error  messages 
Boundarylossmod_passive.f90  -  module  with  array  declarations 

4.1.1  Input  files 

4.1.1 .1  Runinput_v4.inp 

This  is  the  same  file  that  BellhopDRDC_ray_TL_v4  uses,  described  in  section  2.1.1.  The 
boundary  loss  algorithm  uses  the  run  title,  frequency,  and  wind  speed  from  this  file.  For  this 
example,  the  frequency  was  increased  to  3000Hz  and  the  wind  speed  was  set  to  20  kts. 

4.1 .1.2  Speed.inp 

This  is  the  same  file  that  BellhopDRDC_ray_TL_v4  uses,  described  in  section  2.1.2.  The 
boundary  loss  algorithm  uses  the  sound  speed  from  the  surface  and  the  bottom  of  the  first  profile 
in  this  file. 

4.1 .1.3  Bottomloss.inp 

This  is  the  same  file  that  BellhopDRDC_ray_TL_v4  uses,  described  in  section  2.1.3.  The 
boundary  loss  algorithm  uses  all  lines  in  this  file.  For  this  example,  the  acoustic  bottom  was 
described  as  a  simple  half-space  with  speed  =1453  m/s,  density  =  1.41  g/cc  and  attenuation  = 
0.038  dB/mkHz. 
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4.1.2 


Output  files 


4.1 .2.1  Botloss_passive.txt 

This  file  contains  the  bottom  loss  and  grazing  angles  for  each  bottom  description  in  the  input  file. 
The  structure  of  this  file  begins  with  the  run  title,  the  frequency,  the  number  of  range  dependent 
bottom’s  being  specified,  and  a  string  showing  the  type  of  bottom  being  computed,  such  as  the 
two  layer  model,  the  MGS  model  or  a  user  input  table.  Next,  the  number  of  grazing  angles  is 
listed,  followed  by  an  array  of  grazing  angles  then  an  array  of  dB  losses.  These  last  three  are 
repeated  for  each  bottom  specified.  An  example  of  this  file  is  shown  in  Figure  28  generated  using 
a  single  half-space.  An  example  of  the  plot  of  these  data  is  shown  in  Figure  29. 


BELLHOP-  Emerald  basin  toward  Sambro  Bank 


3000.000 

1 


two-layer  fluid  reflection  bottom  loss 
91 

0  12  3  4 

6  7  8  9  10 

12  13  14  15  16 


5 

11 

17 


90 

1.1498389E-I 

9.586963 

21.95951 

31.09075 

22.01334 

19.39529 


1.870405 

11.66344 

25.81959 

27.73483 

21.28639 

19.07889 


3.750788 

13.85982 

31.36469 

25.58814 

20.68758 

18.80259 


5.653327 

16.23368 

42.31104 

24.06392 

20.18612 

18.55959 


7.592482 

18.87825 

37.66504 

22.91468 

19.76052 

18.34455 


Figure  28.  Portion  of  a  Botloss _passive.txt  file. 
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Figure  29.  Plot  ofbotloss _passive.txt  showing  the  single  half-space  region  in  bottomloss.inp  for 
the  acoustic  bottom  descriptions. 


4.1. 2.2  Surfloss_passive.txt 

This  file  contains  the  surface  loss  and  grazing  angles  from  all  available  loss  algorithms.  The 
structure  of  this  file  begins  with  the  run  title,  the  frequency,  a  string  showing  the  surface  loss 
algorithm,  the  number  of  points  in  the  table,  then  an  array  of  grazing  angles  and  an  array  of  losses 
in  dB.  This  format  is  repeated  for  the  second  surface  loss  algorithm.  An  example  of  this  file  is 
shown  in  figure  Figure  30.  A  plot  of  these  data  is  shown  in  Figure  3 1 . 
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Figure  30.  Example  portion  of  surfloss _passive.txt. 


Surface  loss  models 


Figure  31.  Plot  of  surfloss  _passive.txt  for  both  Beckman  Spezzichino  (B)  and  modified  Eckart  (E) 
generated  for  3  kHz  and  20  kt  wind  speed. 


38 


DRDC  Atlantic  CR  201 0-1 34 


4.1.3  Plot  routine 


An  IDL  plot  routine  has  been  prepared  to  provide  a  simple  graphic  representation  of  the 
boundaryloss_passive  outputs.  This  should  be  freely  altered  to  suit  the  users’  data  and  output 
requirements. 

Read  loss  plot  boundaryloss  passive.pro  Routine  to  read  the  output  in 
Botloss_passive.txt  and  Surfloss_passive.txt  and  produce  two  loss  figures. 

4.2  Boundaryloss_active 

Boundaryloss_active.f90  -  main  program 

READINPUT_ACTIVE  -  reads  active  input  files 
READBOTLOSS_S  -  reads  bottomloss.inp 

TMP_SSP  -  function  to  convert  temperature  to  speed  using  Leroy’s  equation 
twolayerRefl  -  computes  reflection  coefficient  from  two  fluid  layers  of  sediment 
CRCIS  -  function  to  convert  real  wave  speed  and  attenuation  to  a  single  complex 
wave  speed 

ERROUT  -  outputs  error  messages 
datamod_active_v4  .f90-  same  as  file  in  BellhopDRDC_active_v4 
refcomod_active_v4.f90  -  same  as  file  in  BellhopDRDC_active_v4 
scatstrength_v4.f90  -  same  as  file  in  BellhopDRDC_active_v4 
SEmod_active  .f90-  same  as  file  in  BellhopDRDC_active_v4 
Sspmod_active_v4.f90  -  same  as  file  in  BellhopDRDC_active_v4 

4.2.1  Input  files 

4.2.1. 1  Active_general.inp 

This  is  the  same  file  that  BellhopDRDC_active_v4  uses,  described  in  section  3.1.1.  The  boundary 
loss  algorithm  uses  the  run  title,  number  of  sensors,  frequency,  wind  speed,  and  ‘EC’  parameters 
p.  n,  and  a  from  this  file. 

4.2. 1.2  Radial_ssp.inp 

This  is  the  same  file  that  BellhopDRDC_active_v4  uses,  described  in  section  3.1.2.  The  boundary 
loss  algorithm  uses  the  sound  speed  at  the  surface  in  the  first  profile  of  the  first  sensor  and  radial 
for  the  Eckart  surface  loss  calculation.  It  also  uses  the  sound  speed  of  the  bottom  in  the  first 
profile  at  each  sensor  and  each  radial  for  the  acoustic  bottom  loss  calculations. 

4.2. 1.3  Radial_bottomloss.inp 

This  is  the  same  file  that  BellhopDRDC_active_v4  uses,  described  in  section  3.1.3.  The  boundary 
loss  algorithm  uses  all  of  this  file. 
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4.2. 1.4  Radial_bathy.inp 

This  is  the  same  file  that  BellhopDRDC_active_v4  uses,  described  in  section  3.1.4.  The  boundary 
loss  algorithm  uses  the  number  of  radials  from  this  file. 

4.2.2  Output  files 
4.2. 2.1  Botloss_active.txt 

This  file  contains  the  bottom  loss  and  grazing  angles  for  each  bottom  description  in  the  input  file. 
The  structure  of  this  file  begins  with  the  run  title,  the  frequency,  the  number  of  bearings  and 
assets,  a  string  showing  the  type  of  bottom  being  computed,  such  as  the  two  layer  model,  the 
MGS  model  or  a  user  input  table,  and  the  number  of  range  dependent  bottom’s  being  specified. 
Next  the  number  of  grazing  angles,  index  of  range  dependence,  bearing  and  asset  are  listed, 
followed  by  an  array  of  grazing  angles  then  an  array  of  dB  losses.  These  last  three  are  repeated 
for  each  bottom  specified  along  the  range,  then  the  output  loops  over  the  number  of  bearings  for 
each  asset.  An  example  of  this  file  is  shown  in  Figure  32.  An  example  of  the  plot  of  these  data 
using  the  IDL  program  read_loss_plot_boundaryloss_active.pro  is  shown  in  Figure  33  for 
bearing  #4  and  asset#  1.  The  different  curves  are  labelled  according  to  their  order  in  the 
radial_bottomloss.inp  file. 
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two— layer  fluid  reflection  bottom  loss 


4.2. 2. 2  Surfloss_active.txt 

This  file  contains  the  surface  loss  and  grazing  angles  from  the  two  available  loss  algorithms 
Eckart  (‘E’)  and  Beckman  Spezzichino  (‘B’).  The  structure  of  this  file  begins  with  the  run  title, 
the  frequency  and  wind  speed,  a  string  showing  the  surface  loss  algorithm,  the  number  of  points 
in  the  table,  then  an  array  of  grazing  angles  and  an  array  of  losses  in  dB.  This  format  is  repeated 
for  the  second  surface  loss  algorithm.  An  example  of  this  file  is  shown  in  Figure  34.  A  plot  of 
these  data  using  the  IDL  program  read_loss_plot_boundaryloss_active.pro  is  shown  in  Figure  35. 
This  output  plot  can  be  compared  with  the  passive  case  output  in  Figure  31,  which  used  a 
frequency  of  3  kHz  and  a  wind  speed  of  20  kts. 
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Figure  34.  Example  portion  of  Surfloss_active.txt. 


Surface  loss  models 


Figure  35.  Plot  of  surfloss._active.txt  for  both  Beckman  Spezzichino  (B)  and  modified  Eckart  (E). 
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4.2. 2. 3  Surfscat.txt 


This  file  contains  the  surface  monostatic  backscattering  strength  and  grazing  angles  from  the  two 
available  scattering  strength  algorithms,  Ogden-Erskine  (‘OE’)  and  Chapman-Harris  (‘CH’).  The 
structure  of  this  file  begins  with  the  run  title,  the  frequency  and  wind  speed,  a  string  showing  the 
surface  scattering  strength  algorithm,  the  number  of  points  in  the  table,  then  an  array  of  grazing 
angles  and  an  array  of  losses  in  dB.  This  format  is  repeated  for  the  second  surface  loss  algorithm. 
An  example  of  this  file  is  shown  in  Figure  36.  A  plot  of  these  data  using  the  IDL  program 
read_loss_plot_boundaryloss_active.pro  is  shown  in  Figure  37. 
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Figure  36.  Example  portion  of  Surfscat.txt. 
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Surface  Scattering  Strength  models 


Figure  37.  Plot  of  surfscat.txt  showing  surface  scattering  strength  models  Ogden-Erskine  (OE) 
and  Chapman- Harris  (CH). 


4.2. 2.4  Botscat.txt 

This  file  contains  the  bottom  monostatic  backscattering  strength  and  grazing  angles  from  the  three 
available  loss  algorithms,  Ellis-Crowe  (‘EC’),  Lambert’s  Rule  (‘LB’),  and  the  Omni  Rule  (‘OM’). 
The  structure  of  this  file  begins  with  the  run  title,  the  frequency,  a  string  showing  the  bottom 
scattering  strength  algorithm,  the  number  of  points  in  the  table,  then  an  array  of  grazing  angles 
and  an  array  of  losses  in  dB.  This  format  is  repeated  for  the  other  two  bottom  scattering  strength 
algorithms.  An  example  of  this  file  is  shown  in  Figure  38.  A  plot  of  these  data  using  the  IDL 
program  read_loss_plot_boundaryloss_active.pro  is  shown  in  Figure  39. 
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Test  run  with  1  sensor,  4  bearings 
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Figure  38.  Example  portion  of  botscat.txt. 
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Figure  39.  Plot  of  botscat.txt  showing  Ellis-Crowe  (EC),  Lambert  (LB)  and  omni  (OM)  bottom 
scattering  strength  models. 


4.2.3  Plot  routine 

An  IDL  plot  routine  has  been  prepared  to  provide  a  simple  graphic  representation  of  the 
boundaryloss_active  outputs.  This  should  be  freely  altered  to  suit  the  users’  data  and  output 
requirements. 

Read  loss  plot  boundaryloss  active.pro  -  Routine  to  read  the  output  in  Botloss_active.txt, 
Surfloss_active.txt,  surfscat.txt  and  botscat.txt  and  produce  four  loss  figures. 
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